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ABSTRACT 



Initial steps are reported in the development of a method for deter- 
mining the instantaneous velocity vector at the exit of a compressor rotor 
using two Kulite pressure probes with synchronized sampling. A method 
to determine the yaw angle of the flow at any point in the reference frame 
of the rotor using a single cylindrical Kulite impact probe and an equivalent 
pneumatic probe is described. The method was verified in the steady flow 
of a free jet and applied in the unsteady flow field of the compressor. On- 
line calibration of the Kulite probe with reference to the pneumatic probe 
was essential to the success of the method. Results for the distribution 
of the yaw angle peripherally across the exit of chosen pairs of rotor blade 
passages, are given. The analytical approximation of impact probe 
characteristics, on which the method is based, and its application with 
experimental data are described. 
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I. INTRODUCTION 



The work reported here is part of an on-going effort at the Naval 
Postgraduate School to determine the aerodynamic characteristics and 
performance of a single stage, axial transonic compressor (Fig. 1) using 
real time instrumentation. 

In work already reported a multiple sensor pneumatic probe (Dodge 
probe) was developed to determine the average velocity at the rotor exit 
[Ref. l]. The design of the probe was based on a knowledge of the charac- 
teristics of simple impact probes, and experience which had been gained 
in representing analytically the characteristics of multiple sensor probes 
[Refs. 2 & 3 ]. A new method of representing the characteristics of mul- 
tiple sensor probes was developed and the time averaged flow from the 
rotor was measured [Ref. 4]. In order to obtain information on the flow 
within the rotor itself, real-time measurements were required. A syn- 
chronized sampling system was therefore developed which allowed program- 
mable digital data acquisition from fixed instrumentation. [Ref. 5] Using 
this electronic "pacer, " sampling and analog to digital conversion of trans- 
ducer signals was controlled to be at any of 128 positions between any 
selected adjacent pair of the 18 blades of the rotor, independent of the 
speed. 

A method of velocity measurement was proposed which used two Kulite 
probes with synchronized sampling. The technique would establish the 
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instantaneous and some time dependent properties of the velocity field in 
the frame of the compressor rotor. [Ref. 6] In addition it was proposed 
to relate the results obtained with this method to measurements made 
using hot wire probes and using laser velocimetry. The primary objec- 
tive of the overall research effort was to understand and to interpret the 
phenomena occurring at different machine conditions. [Ref. 6] 

The purpose of the present work was to begin the development of the 
proposed probe system to measure "instantaneous” velocity at the exit 
of the transonic compressor rotor. The ideas behind the two-probe tech- 
nique were derived from the previous work on multiple sensor probes 
[Ref. 1,2,3&4] and the application of synchronized sampling. Concept- 
ually the two probes can be considered to supply measurements which 
correspond to those of the sensors of a multiple sensor probe. The 
synchronized sampling allows measurements from two physically separated 
probes to be taken (at different times) at the same point in the reference 
frame of the rotor. 

The initial step in the proposed method was to determine the yaw 
angle of the flow at any required point in the rotor frame. In the present 
work a method for determining the flow yaw angle at a point in the rotor 
frame using a single Kulite probe was developed and successfully tested 
in the transonic compressor. 

Section II describes the methods used to obtain velocity from pressure 
measurements taken with a multiple sensor probe, and introduces the 
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proposed velocity measurement method and the method to determine flow 
yaw angle. In Section III the Kulite and pneumatic probes and the data 
acquisition system used in the current work are described. A technique 
for calibrating the Kulite probe on-line to an equivalent pneumatic probe 
was developed and the method for yaw angle measurement was verified in 
experiments using a steady, calibrated free jet. The tests and results 
are presented in Section IV. Section V describes the transonic compres- 
sor test conducted to determine the flow yaw angle at selected points in 
the rotor frame and discusses the results. Conclusions of the study are 
given in Section VI. 

The use of an impact probe to measure yaw angle required that the 
characteristics of such a probe be known. Measurements of the character- 
istics of an impact probe were made in a free jet, and an analytical expres- 
sion was found to represent the results. The analytical representation 
was fundamental to the successful measurement of the yaw angle, which 
followed. The measurements and analysis of the results from tests of a 
simple impact probe are given in Appendix A. 

The derivation of a numerical procedure to approximate cylindrical 
probe characteristics with the empirical expression is given in Appendix B. 
Als© described is the method used to derive yaw angle from probe data 
using the same empirical expression. Finally, the empirical expression 
was used to derive the calibration characteristics of the Dodge probe 
which was constructed using tubing similar to that used in the construction 
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of the test impact probe. The derivation is given in Appendix D. Ap- 
pendix C describes the computer programs used in acquiring data. 
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II. APPROACH 



A. USE OF MULTIPLE SENSOR PROBES TO MEASURE VELOCITY 
The United Sensor five hole probe and the probe constructed by 
F. J. Dodge [Ref. l], shown in Figure 2 can be used to measure air 
average velocity when the outputs of the multiple sensors oriented at 
various angles to the probe axis are calibrated over a range of Mach 
numbers and pitch angles. Yaw angle is first determined directly by 
balancing the outputs of the two sensors which are set at equal and op- 
posite angles to the plane of the impact sensor and the probe shaft. Pitch 
angle (0) and non-dimensional velocity magnitude (x) are determined from 
readings taken after the probe has been rotated to balance the pressures 
at sensors Pj^ and P2 (Fig. 2 ), and the yaw angle has been noted from a 
vernier scale. 

The first method for calibrating and applying multiple sensor 

probes uses the definitions: 

^ (x, 0 ) = (Pj - ^ 23 ^ ^ ^1 (United Sensor and Dodge) 

7 = 7 (X, 0 ) = (Pj - P4) / (Pj - P23) (Dodge) 
or 

If - If = (P^ - P5) / (Pj - P23) (United Sensor) 



^Footnote: The non-dimensional velocity x is defined as x = V/V^, where 

V = velocity magnitude and = "limiting velocity" = 2 CpT 
for a perfect gas, where Cp = specific heat at constant ^ 

pressure and Tj- = stagnation temperature. 
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In the calibration, data for j2 and 'J are measured in a controlled 
flow over the range of pitch angles and Mach numbers expected in the un- 
known flow to be measured. The data is then approximated to yield fifth 
order polynomial expressions for, '~jf = 7 ^( 0 , x) and x = x( J2 , 0 ). In 

application, probe measurements in an unknown flow are reduced to the 
velocity vector by solving the polynominals simultaneously. This method 
is described in Refs. 1, 2, and 3. 

A second calibration method [Ref. 4] has been developed in which 
and are defined as in the first method. Calibration data is obtained in 
a known flow and reduced to fifth order polynomial expressions for: 



0 = 0 ( y , x) 
fi = {<i, x) • U (x) 

where U {x) - { ^ I ^ -1) x. (1- x^) 

From these approximations, which include the primary dependence 
on velocity explicitly, readings of and If in an unknown flow are re- 
duced to velocity magnitude and pitch angle as in method 2. 

A third method uses an analytical expression for the characteristic 
behavior of a single cylindrical impact tube in yaw or pitch at different 
Mach numbers. The overall calibration for a multiple sensor probe 
constructed using similar tubing then can be derived analytically as shown 
in Appendix D. 
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B. MEASUREMENT OF VELOCITY USING 2 KULITE SEl^ICONDUCTOR 
PROBES WITH SYNCHRONIZED SAMPLING 

Construction of sufficiently small multiple sensor probes using Kulite 
semiconductor sensors was not considered to be practical due to the small 
probe size required and the low probability that a probe could be construc- 
ted without the failure of at least one sensor. Moreover, such a probe 
could not be used to resolve peripheral variations in the yaw angle down- 
stream of the transonic rotor. 

What was proposed in Ref. 6 was to use two separate impact probes 
to perform the tasks of the separate sensors in the multiple sensor probe. 

The two probes are located at different peripheral but identical axial 
and radial stations between the rotor and stator blade rows. The two pro- 
bes are shown in Fig. 3. The probes are sampled under pacer control 
such that, at the time the two readings are digitized, the probes are at 
the same point in the reference frame of the rotor. [Ref. 6 ] A computer 
peripheral device which is required to control the data acquisition (pacer) 

was developed by West. [Ref. 5] 
o 

The 90 or type A probe shown in Figure 3 can be rotated in yaw 

about the tip to angles corresponding to the sensors for Pj, ^3 

o 

in the Dodge probe. The 55 or type B probe shown in Figure 3 must be 

« 

set at the balanced or zero yaw angle (determined using the type A probe) 
to measure the pressure corresponding to in the Dodge probe. To 
determine velocity, one of the three methods described in the preceding 
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section must be used to represent and apply the calibration. The two 



probes must be calibrated together in a steady uniform flow as if they 
were the sensors of a multiple sensor probe. 

The initial step in applying the two probe system to measure velocity 
is to determine the yaw angle of the flow at the required point in the rotor 
frame. Only then can the two probes be set to angles at which the cali- 
bration was established. In the present work, a method for determining the 
yaw angle at each point in the rotor frame was developed using only the 
type A probe. 

Since the method makes use of the characteristics of cylindrical impact 
probes, a preliminary investigation was made of the characteristics of 
pneumatic cylindrical impact probes. The resvilts are given in Appendix A. 
It was found that the response of the probe could be represented as a pres- 
sure coefficient Cp, defined as: 




N 



where 



A, B and N are constants for the particular probe. 



M = Mach number 



p = static pressure 



= ratio of specific heats 



= probe pressure 



Pj. = impact pressure 
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= angle of the flow to a reference scale 
angle of the axis of the probe to the reference scale 
When the flow angle to the probe axis was a combination of pitch angle, 
0, and yaw angle, OC , then was calculated using the geometrical 
relationship: 

= cos “^(cos 0 cosOC ) 

These results provided the basis for the development of a technique 
to derive yaw angle from Kulite probe measurements: the development is 
described in Section IV, following a description of the probes and data 
acquisition method given in the next section. 
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m. INSTRUMENTATION 



A. PROBE DESCRIPTION 

Figure 3 is drawing of the XB-062-25 semiconductor probes manufac- 
tured by Kulite for the tests conducted in this study. The sensors in each 
probe tip were standard CQ-052 series Kulite ultra miniature pressure 
sensors. Kulite Type "B” screens were installed on the probes to protect 
the sensors from particle impingement in high velocity flows. Figure 4 
shows the probe tip dimensions and internal design. 

Kulite -equivalent pneumatic probes were constructed at Naval Post- 
graduate School that were similar (to within 0. 002 inches) to the external 
geometry of the Kulite probes shown in Figures 3 and 4. Initial impact 
pressure measurements taken with the Type A Kulite and pneumatic 
equivalent probes were consistently 10% below actual total pressure in 
the flow. It was immediately recognized that flow stagnation was occurring 
in the center of the probe face and not at the holes in the screen. There- 
fore, a shroud of 0.083 inch outside diameter, stainless steel tubing was 
added as shown in Figure 4. Subsequent tests indicated good agreement of 
impact pressure measurements with actual total pressure with no notice- 
able loss in the frequency response of the Kulite probe. 

B. DATA ACQUISITION SYSTEM 

Data acquisition from the Kulite probes and conventional strain gauge 
transducers was controlled by a Hewlet-Packar d 21MX computer. A pacer 
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developed by West [Ref. 5] was used for "synchronized" sampling behind 
the transonic compressor rotor as shown in Figure 5. Synchronized 
sampling allows the properties of the flow to be determined at a point 
which is fixed with respect to the rotor. The pacer uses conditioned 
1 per blade and 1 per rev signals to generate a trigger pulse at a selected 
location in the rotor frame. There are 256 equally spaced selectable 
locations in each rotor blade pair. Steady and nonsteady flow properties 
can be determined by sampling at a discrete location as many times as 
is necessary. Figure 6 shows the complete data acquisition system used 
in the present study. Appendix C is a discussion of the software used 
with the 21 MX computer. 
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IV. PRELIMINARY MEASUREMENTS IN A FREE JET 



A. INTRODUCTION 

In order to use semiconductor probes for quantitative pressure meas- 
urements in an environment in which the temperature is unknown, a means 
of on line calibration must be used. The output voltage of the Kulite probe 
is linear with respect to the pressure difference across the sensing 
diaphragm. A typical Kulite calibration is shown in Figure 8. However, 
large shifts in output voltage due to temperature changes have been re- 
ported and were confirmed by Paige [Ref. 7]. For this reason, Kulite 
semiconductors have not been used for absolute pressure measurements 
in turbo machinery but have been used to measure the fluctuations in un- 
steady or periodic flows. In calibration tests conducted in a steady air 
flow early in the present work, these shifts of level of output were ob- 
served. Shifts in the slope of the output were also observed when the 
Kulite probe was rotated to a new flow angle. Therefore, a technique to 

calibrate the Kulite probe against a geometrically identical pneumatic 

\ 

probe, similarly oriented in the flow, using the computer on line, was 
developed. 

The goal of the experiments described in this section was to develop 
the on-line calibration procedure and to verify, in steady flow, a method 
for measuring yaw angle in the compressor. A four inch free jet was 



20 



used to provide a known uniform flow field for the tests. The test appara- 
tus is described in reference 8. The Kulite type A and the equivalent 
pneumatic type A probes were inserted in the flow to be two inches apart 
on opposite sides of the jet axis. The mounting apparatus allowed both 
probes to be pitched and yawed simultaneously and to be set to similar 
angles with respect to the flow direction. A Prandtl probe was used to 
monitor the jet velocity. 

B. KULITE-EQUIVALENT PNEUMATIC PROBE TEST 

Tests in the free jet, similar to those described in Appendix A, were 
conducted to establish the characteristics of the type A equivalent pneu- 
matic probe with respect to yaw angle and Mach number. The range of 
Mach numbers surveyed was . 4 to .6. 

The results are shown in Figure 7. The pressure coefficient vs. 
yaw angle characteristics for this probe and the variation with Mach 
number were found to be qualitatively similar to those of the cylindrical 
impact probe reported in Appendix A. The coefficients A, B, N, and 
were successfully calcvilated as described in Appendix B. It is concluded 
therefore that the technique for calculating the zero yaw angle described 
in Appendix B could also be used with the type A probe geometry. 

C. KULITE PROBE TESTS 

1 . Averaging Techniques 

The first measurements from the Kulite probe in a free jet veri- 
fied the high level of turbulence and unsteadiness which were known to be 
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present from earlier measurements. In order to investigate the charac- 
teristics of the Kulite probe with respect to yaw angle and Mach number, 
some method of averaging over a number of data samples was necessary. 

The purpose of the first test therefore was to determine the effect 
of sample number and interval on the average measurement obtained from 
the Kulite probe. 

At one test condition, three ten minute tests were conducted in 
which the techniques for sampling the Kulite probe were varied. During 
the tests, the Kulite equivalent pneumatic probe pressure was read at 
close intervals using a water column U-tube manometer referenced to 
atmosphereic pressure. The Kulite reference pressure was sensed by 
a conventional strain gauge transducer connected to and calibrated on 
one input channel of the A/D converter. The reference pressure was 
obtained at one minute intervals by recording the average 100 samples 
taken at 10 millisecond intervals. 

In each of the three tests, a data point for the Kulite probe 
pressure was recorded at 1 minute intervals for a period of 10 minutes. 
Each data point consisted of taking the average voltage from a number 
of ensembles of samples taken at 10 microsecond intervals. The data 
for the methods used in the three tests is shown in the following table: 
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METHOD 



SAMPLES PER 
ENSEMBLE 



ENSEMBLE 
# OF INTERVAL 

ENSEMBLES (SEC) 



TOTAL TIME 
(SEC. ) 



1 


100 


1 


- 


1 


2 


1000 


1 


- 


10 


3 


100 


10 


2 
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For each data point, the pressure was calculated using the 
average voltage from the Kulite probe and the recorded reference pres- 
sure as shown in the following section, C.2. The results are shown in 
Figure 9. 

It can be seen that methods two and three gave acceptable accuracy 
whereas method one did not. The failure of method one was the result of 
high frequency turbulence in the free jet. Seen on the oscilloscope, the 
dominant fluctuations in the Kulite signal had a period of about two milli- 
seconds, so that a sample duration of one millisecond could not give a 
correct time average. Method two was used in subsequent tests in the 
free jet. 

2. Calibration Method 

The need to calibrate the Kulite transducer to account for changes 
due to temperature, and for shifts in output slope which occurred when 
the probe was rotated in the flow, was stated in Section IVA. A calibra- 
tion procedure was therefore developed in the free jet tests and subsequen- 
tly applied in the measurement of yaw angle. The Kulite probe calibra- 
tion used the assumption that the equivalent pneumatic probe measured 
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the time average of the pressure on the face of the Kulite probe when the 
two probes were similarly oriented in the flow. The calibration procedure 
was to apply a controlled pressure to the reference tube of the Kulite 
probe, then sample, (as described in Section IV C. 1) the pressure from 
the equivalent pneumatic probe, Pp, the Kulite reference pressure, Pr, 
and the Kulite output voltage. This procedure was repeated for a number 
of calibration points, changing the reference pressure for each point. The 
data was reduced on-line using a linear least squares routine to obtain 
the constant coefficients and xj in the equation, 

Pp - Pr = Xo + xj 1E 

where E was the average amplified output voltage of the Kulite transducer. 

In subsequently applying the calibration, samples of Kulite output 
voltage were reduced on line to values of pressure, Pk, using the calibra- 
tion coefficients x^ and Xj in the equation 
Pk = Xq + Xj^ • IT + Pr 

An experiment was conducted to determine the minimum number 
of reference pressure settings required to calibrate the Kulite probe with 
acceptable accuracy. A calibration test was carried out, as described 
above, using eight different reference pressures. The coefficients x^ and 
x^ were computed first using all eight data points. Then, using the same 
data, new coefficients were computed using various subsets of the eight 
data points. Table IV- 1 shows the subsets used and the deviation of the 
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results from the eight point calibration. The results indicated that only- 
small improvements in the accuracy were obtained by using more than 
two points. Two points, if properly chosen to include the range of ex- 
pected pressures, gave a sufficiently accurate calibration. In subsequent 
experiments a two-point calibration was therefore used with atmospheric 
pressure as one reference pressure and the flow total pressure as the 
other. 

D. YAW ANGLE DETERMINATION 

The purpose of this experiment was to verify, in steady flow, a method 
of measuring the zero yaw angle behind the rotor of the transonic compres- 
sor using the type A Kulite probe. Measurements were made in the four 

inch free jet at Mach numbers of 0.406 and 0. 587. At each Mach number 

o o o 

both probes were set at pitch angles of -10 , 0 , and +10 . At each pitch 
angle both probes were set at yaw angles of +60 , +45°, +30°, +15°, and 
0° to the flow. After each change of angle the Kulite probe was calibrated 
as described above. Water column readings of the Prandtl probe and the 
Type A pneumatic probe were recorded, the Type A pneumatic pressure 
was sampled and the Type A Kulite average voltage was sampled and 
reduced on line to values of pressure. 

In the off-line analysis, zero yaw angle was calculated from the data 
using the least squares method described in Appendix B. The number of 
points needed to determine the zero yaw angle was investigated by varying 
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the data points included in the least squares calculation. The results 
are given in Table IV-2 for 3, 5, 7 and 9 points. 

It was seen that the zero yaw angle was determined to within 0. 5 
degrees for pitch angles of -10, 0 and +10 degrees. It was concluded that 
the Kulite probe, when carefully calibrated to an equivalent pneumatic 
probe, can be used to measure yaw angles with reasonable accuracy. 

The need to calibrate the Kulite probe at each orientation was con- 
firmed by an inspection of the calibration coefficients determined in the 
experiment. The coefficients are shown in Table IV-3. 

The variation in the coefficients appeared to be sufficiently random 
in nature as to preclude any attempt to describe the change as a function 
of Mach number and flow angle. The variations were not understood, 
however, the technique of recalibrating at each angle was one which could 
also be applied in the compressor. It should be noted, however, that in 
the periodic flow behind a compressor rotor, application of the same 
calibration and measurement techniques will force agreement between 
the absolute values of the pneumatic time average and Kulite time average 
pressures. 



26 



(i) DATA 



POINT P-Pr water) AMPLIFIED KULITE OUTPUT 

(millivolts) 



1 


1 

• 

00 




-68. 9415 


2 


-27. 5 




-46. 193 


3 


-12. 5 




-25. 2698 


4 


- 0.2 




- 6.05109 


5 


24. 2 




29. 0726 


6 


36. 0 




46. 186 


7 


47. 2 




62.0352 


8 


57. 8 




77. 7253 


Impact pressure at the probe face, p^= 


37.9 ins. water gauge. 


Least squares fit to linear 


calibration 


given by, P -Pj.=Xq+Xj^E : 


x^=4. 39565 ins. 


water, 3Tj^ = 


=686. 358 ins 


. water /volt. 


(ii) Least Squares fit to linear calibration, P-P =x +x E, using 


subsets of points. 






Points used for 


Subset 




(Xj-Xj)/Xj (%) 


1.3, 5,7 




0. 2 


0. 16 


00 




0. 1 


0. 01 


2,4,7 




-0. 1 


-0. 60 


3,4,6 




0. 0 


0. 70 


1,3,5 




0. 1 


0. 44 


3,5,7 




0. 3 


0. 50 


1,8 




0. 2 


0. 07 


2,7 




0. 0 


-0. 51 


3,6 




0.3 


1. 14 


2,5 




-0. 1 


-0. 36 


1.5 




0. 0 


-0. 10 


2, 6 




-0. 1 


0.39 


3,7 




0. 4 


1 

o 

o 


4,6 




-0. 2 


-0. 11 



Table IV- 1 Results of Kulite probe calibration in. steady flow 



im 

1 




Pj^=Kulite pressure calculated using linear calibration 

Pj^^=Kulite equivalent pneumatic pressure measured with conventional 
transducer 

(P, ) =Kulite equivalent pneumatic pressure measured with water 

ke m ^ ^ 

column manometer 

Entries in the table are the zero yaw angles ( q), in degrees, calculated 

from the above measurements using subsets of the data recorded at a total 
of nine angles. 

The subsets are as follows; 



Subset # Angles included 





1 




0°,+15°,+30°, 


+45°, + 60 ° 






2 




0°,+30°,+45°, 


+60° 






3 




0°,+45°,+60° 








4 




0 °,+ 60 ° 






Subset 




No. of points Q(_ ^from Pj^ 


(X from P, 
^o ke 


a, from 








(M=0.406, 0=0°) 






1 




9 


.70 


. 80 


. 81 


2 




7 


.75 


. 80 


. 81 


3 




5 


. 82 


. 86 


. 86 


4 




3 


1.31 


1. 29 


1.33 








(M=0. 4O6,0=+1O°) 






1 




9 


1.38 


1. 10 


1.08 


2 




7 


1.40 


1.11 


1.09 


3 




5 


1. 48 


1. 15 


1. 14 


4 




3 


2.09 


1.61 


1. 60 








(M=O.4O6,0 = -1O°) 






1 




9 


.46 


. 60 


.60 


2 




7 


.49 


.62 


.60 


3 




5 


.60 


. 65 


.62 


4 




3 


.78 


.98 


00 

• 



Table IV-2 Calculated zero yaw angles from test described 

in Section IV -D. 
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Subset 



No. of points 



oc 



from P, 
o k 



a^from 



(P^^) 



(M=0. 587, 0=0°) 



1 


9 


.86 


.85 


.85 


2 


7 


.86 


.85 


. 86 


3 


5 


.95 


.93 


.92 


4 


3 


1.33 


1.47 


1.44 






(M=0.587, 0 = 10°) 






1 


9 


.91 


.79 


. 83 


2 


7 


.93 


.81 


.84 


3 


5 


.95 


.84 


. 88 


4 


3 


1.38 


1. 17 


1.25 






(M=0. 587, 0=-lO°) 






1 


9 


. 85 


.77 


. 73 


2 


7 


. 85 


.78 


.74 


3 


5 


.89 


.82 


.78 


4 


3 


1.35 


1. 28 


1. 23 



Calculated zero yaw angles from test described 
in Section IV-D. 



lV-2 (cont) 




I 





0=0° 




0=- 


0 

10 


0 = 


+10° 


robe Yaw 


X 


^1 


X 


^1 


X 


^1 


ngle (Deg) 


o 


0 

(M=0. 406) 


0 


1 


0 


-3.9 


696 


-3.9 


686 


-5. 8 


685 


15 


-5.5 


721 


-4.2 


694 


-5. 9 


691 


-15 


-3.7 


686 


-5.2 


691 


-6.3 


694 


30 


-2. 6 


683 


-2. 8 


693 


-5. 1 


690 


-30 


-5.0 


689 


-5.3 


692 


-7.4 


691 


45 


-0. 5 


685 


-0. 2 


689 


-2.3 


694 


-45 


-7.0 


686 


-7.4 


698 


-9.4 


677 


60 


+4. 5 


715 


+3. 6 


681 


+3. 1 


679 


-60 


-9.3 


711 


-9.6 

(M = 0. 587) 


693 


-12.4 


671 


0 


-4.8 


689 


-5. 0 


692 


-4.5 


685 


15 


-4.0 


690 


-5.0 


694 


-5.3 


684 


-15 


-6. 5 


691 


-6. 1 


689 


-6. 1 


687 


30 


-3.5 


692 


-1.3 


697 


-3.8 


681 


-30 


-9. 1 


692 


-8.5 


690 


-10. 0 


691 


45 


-2. 0 


735 


+5. 5 


684 


+0. 9 


692 


-45 


-11. 1 


688 


-12. 6 


685 


-13.7 


688 


60 


+16. 1 


673 


+16.7 


684 


+16,3 


672 


-60 


1 

o 

• 

o 


697 


-20.4 


680 


-19.6 


659 



X in inches water 
o 

in inches water/volt 



e. 



Table IV-3 Variation of Kulite calibration coefficients in 
yaw angle test on the free jet. 
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V. FLOW ANGLE MEASUREMENTS IN A TRANSONIC COMPRESSOR 



A. TRANSONIC COMPRESSOR 

The transonic compressor is shown in Figure 1. It is driven by an 
air turbine drive unit capable of supplying 450 horsepower at 30, 000 RPM. 
When operating at the design point, the relative Mach number at the rotor 
blade tip is 1. 5. The flow rate is controlled by an electrohydraulic rotat- 
ing throttle plate located at the inlet duct, which also contains a filter and 
flow measuring nozzle. An Allis -Chalmers multistage axial compressor 
supplies the turbine drive air. A complete description of the test facilities 
is given in reference 8. 

B. PROBES AND INSTRUMENTATION 

A Kulite type A probe, a pneumatic type A probe (Fig. 3), and the 

Dodge probe (Fig. 2) were inserted into the compressor at an axial distance 

of 0.65 inches (.31 chord lengths) behind the rotor trailing edge. A similar 

probe installation is shown in Figure 10. The radial displacement of 

each probe tip from the case wall was 0. 9 inches. The three probes were 

o 

separated peripherally at 45 intervals. 

The pneumatic connections to the Dodge probe PI sensor and to the 
pneumatic type A probe are shown in Figure 11. The pneumatic pressures 
were sensed by conventional differential strain gauge transducers referenced 
to atmosphere. 
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The two strain gauge transducers and the Kulite type A probe trans- 
ducer were connected through conditioning circuits to separate input 
channels of the high speed data system shown in Figure 6, The two 
transducers were calibrated using a water column manometer as a 
standard. The calibration coefficients obtained were used in the on-line 
program to reduce the transducer voltage outputs to pressure in inches of 
water gauge. A control program was written and used with the program 
"DATACQ” to acquire data during the tests. These programs are des- 
cribed and listed in Appendix C. On-line calibration of the Kulite trans- 
ducer was an integral part of the experimental procedure. 

C. TEST PROCEDURE 

The compressor was stabilized at the desired operating condition and 

steady state performance data were recorded. The Dodge probe was 

rotated to balance the pressures at sensors P2 and P3 and the yaw angle, 

QC , was recorded. This angle, 0^ , was taken as a reference angle 

for the Type A probes. The pneumatic and Kulite Type A probes were 

rotated together to the following angle settings, in turn: CX»CXi ^0° , 

~ o — o 

0(i. OCi. each setting, the following procedure was carried out: 

1. The Kulite probe output was sampled in Pacer "Free Run" mode. 
1681 samples were taken at 10 microsecond intervals, the aver- 
age was calculated and reduced on line to pressure. 

2. The pressure at the PI sensor of the Dodge probe and the pneu- 
matic equivalent pressures were sampled. 
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3. Kulite output was sampled in a survey at 128 locations across 
blade pair 2 using synchronized sampling. Each location was 
sampled 10 times., (on successive revolutions) and the average 
of the 10 samples was reduced online to pressure. 

4. The pressure survey data were displayed on the HP 9862A plot- 
ter and stored in the HP9768 mass memory system. 

5. Steps 3 through 5 were repeated for blade pair 8. 

Data recorded during the test are given in Table V-1. Plots of the 
data obtained by synchronized sampling are shown in Figures 12 and 13. 
The synchronized data was observed to agree qualitatively with oscil- 
loscope traces of the Kulite output voltage. 

D. DATA REDUCTION TO FLOW ANGLE 

The data from seven angles were used to calculate a zero yaw angle 
for each discrete location within the blade pairs, using the method des- 
cribed in Appendix B. The BASIC prdgram "KAW78" used to reduce the 
data is listed in Table V-2. 

E. RESULTS AND DISCUSSION 

The distributions of the zero yaw angle for the Kulite probe (the yaw 
angle of the flow relative to the axial direction) for blade pairs 2 and 8 
are shown in Figures 14 and 15 respectively. 

In both figures, the wakes of the rotor blades are evident, and fluctua- 
tions are seen to be present outside of the wakes. It should be noted that 
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the shape of the observed angle variation through the wake would be the 

result, qualitatively, of the flow velocity being non-uniform in magnitude 

but at a constant angle relative to the rotor. The fluctuations outside 

the wakes are not understood. They are seen to be more definite in blade 

pair 2 than in blade pair 8. Further analysis of the data and more data 

at different test conditions are needed to determine the accuracy of the 

detail which is evident in Figures 14 and 15. 

The numerically averaged yaw angles for blade pairs 2 and 8 were 
o o 

22. 8 and 22. 5 respectively. The pneumatically averaged yaw angle 

o 

obtained from the Dodge probe was 23.2 . The yaw angle calculated using 

o 

time averaged pressures from the equivalent pneumatic probe was 22.2 . 
This agreement supports the magnitudes of the yaw angles shown in 
Figure 14 and Figure 15, for which no other verification is presently 
available. 
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Compgressor Performance Data 

N = 18360 RPM 

Weight Flow = 11.92 Ib/sec 



Ttl = 


76. 1 


F 
















cr 

II 


389. 


7 inches 


water 














Pres! 


sure 


Ratio (t- 


t) = 1. 


167 












Efficiency 


(t-t) = 0 


.92 














Data 


From Probe ^ 


Measurements 












Presi 


sures 


in inches of water 
















Blade 


Probe 


X 

o 


^1 


Pk-Pa 


^ka' 


■Pa 


Pl-Pa 


Point 




Pair 


Angle 




















in. 


in. 








(Dodge 








(deg) 


(H^O) 


(H^O/volt) 






Probe) 


1 




2 


23. 2 


-4.02 


702 


55. 0 


56. 


2 


56. 4 


2 




8 


23.2 


- 


- 


56.9 


- 




- 


3 




2 


53. 2 


-6.44 


689 


50. 2 


51. 


2 


57. 1 


4 




8 


53. 2 


- 


- 


52. 9 


- 




- 


5 




2 


-6.8 


-0.90 


691 


52. 7 


51. 


5 


55. 7 


6 




8 


-6.8 


- 


- 


51.3 


- 




- 


7 




2 


68. 2 


-14. 34 


702 


31.9 


32. 


6 


56. 4 


8 




8 


68. 2 


- 


- 


32.2 


- 




- 


9 




2 


-21. 8 


2. 20 


684 


37. 1 


3 7. 


2 


56. 6 


10 




8 


-21.8 


- 


- 


37.2 


- 




- 


11 




2 


78.2 


-14. 77 


666 


18. 8 


19. 


8 


56.0 


12 




8 


78. 2 


- 


- 


19.7 


- 




- 


13 




2 


-31.8 


9.44 


731 


18.4 


21. 


6 


55. 3 


14 




8 


-31. 8 






19. 1 









TABLE V-1 TRANSONIC COMPRESSOR TEST DATA 
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TABLE V-2. BASIC PROGRAM "KAW78'' 
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TABLE V-2 (continued) 
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710 NEXT I 
720 NEXT J 
725 MAT PRINT 
730 RETURN 



VI. CONCLUSIONS 



In measurements made downstream of the transonic compressor 
rotor, there was found to be good agreement (better than 1%) between the 
yaw angle determined using a Dodge pneumatic probe and the average of 
the peripheral distribution of the yaw angle calculated from a system of 
Kulite probe measurements. It was therefore concluded that the tech- 
nique to determine "instantaneous” yaw angle proposed here and verified 
in a steady free jet, also gave good results in the unsteady compressor 
flow field. The yaw angle variation in the rotor blade wakes was resolved 
and the results contained very little scatter. A weak periodic structure 
was detected across one of the two pairs of blade passages measured 
(Blade pair #2). More analysis of the data is required however, before 
conclusions can be drawn concerning the detail which appears to be present 
in the results. 

The yaw angle measurement was a necessary step in the development 
of a method to determine the velocity field at the rotor exit from stationary 
instrumentation. The next step is to determine the best method for repre- 
senting and applying the calibration of a two probe system. The fluctua- 
tions measured in the zero yaw angle (or turning an§le) behind the trans- 
onic compressor rotor makes less attractive, but does not preclude, a 
calibration method that requires the type B probe to be set at the zero 
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yaw angle for each location sampled. However, an interactive computer - 
probe control system would be needed. Alternatively, a method based on 
the approach developed in Appendix D has potential for success if applied 
to the proposed two probe system. The multiple -sens or interference 
effects found in the Dodge probe, to which the method of Appendix D was 
applied, are not present in the two probe system. 

The development of such a method would require that suitable expres- 
sions be found which would individually characterize the behavior of the 
type A and type B probes with respect to the Mach number, yaw angle, 
and pitch angle of the flow to be measured. The expression developed in 
this study was shown to characterize properly the behavior of cylindrical 
probes with respect to Mach number and yaw angle, but the dependence 
on pitch angle was not obtained explicitly. 

Reasonable agreement was noted (Table V-1) between the magnitudes 
of the time averaged Kulite probe impact pressure and the pneumatically 
averaged impact pressure obtained using the Dodge probe. It was con- 
cluded that a Kulite probe, with on-line calibration to an equivalent pneu- 
matic probe, can be used for quantative real time impact pressure measure- 
ments downstream of a rotor. 
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FIGURE 1. TRANSONIC COMPRESSOR 
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STRAIGHTENER & kULITE 
TRANSDUCER 
PORTS ( 18 ) 
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FIGURE 2. DODGE COMBINATION PROBE PHOTOGRAPH WITH 
THERMOCOUPLE REMOVED 
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FIGURE 3. KULITE PROBES 
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FIGURE 4. PROBE TIP DETAIL 
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FIGURE 7a 

CHARACTERISTICS OF THE KULITE -EQUIVALENT PNEUMATIC PROBE 
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FIGURE 7b 

CHARACTERISTICS OF THE KULITE-EQUIVALENT PNEUMATIC PROBE 
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FIGURE 7c 

CHARACTERISTICS OF THE KULITE-EQUIVALENT PNEUMATIC PROBE 
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FIGURE 8. KULITE PROBE CALIBRATION CURVE 
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FIGURE 9 

EVALUATION OF AVERAGING TECHNIQUES IN THE KULITE PROBE 

CALIBRATION TEST 
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FIGURE 10 

PHOTOGRAPH OF TYPE B PNEUMATIC EQUIVALENT PROBE 
IN COMPRESSOR MOUNTING APPARATUS 
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FIGURE 11. DIAGRAM OF COMPRESSOR TEST INSTALLATION 
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FIGURE 12b 

BLADE PAIR 2 KULITE PRESSURES 
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BLADE PAIR 2 KULITE PRESSURES 
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BLADE PAIR 2 KULITE PRESSURES 
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FIGURE 13a 

BLADE PAIR 8 KULITE PRESSURES 
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FIGURE 13b 

BLADE PAIR 8 KULITE PRESSURES 
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FIGURE 13c 

BLADE PAIR 8 KULITE PRESSURES 
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FIGURE 13d 

BLADE PAIR 8 KULITE PRESSURES 



63 



PRCER 5RMPLE POINT 




PK-Pfl (IN M) 



T 

N 



S 

P4 

rvi 



s 

s 

w 



IS 

m 



s 

in 



s 

r 



s 



s 

IS 



s 

m 



s 

in 



s 

T 



S 

ru 



IS 



FIGURE 13e 

BLADE PAIR 8 KULITE PRESSURES 
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FIGURE 14. BLADE PAIR 2 ZERO YAW ANGLES 
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APPENDIX A 



DETERMINATION OF THE CHARACTERISTICS OF 
CYLINDRICAL IMPACT PROBES 

A. 1 Calibration Tests 

A cylindrical impact probe was constructed to have a tip with the 
same dimensions as the individual sensors in the Dodge probe. The probe 
is shown in Figure Al. Calibration tests of the probe were conducted in a 

seven inch free air jet at five Mach numbers varying from 0, 2 to 0. 5, 

o o o o 

yaw angles from -90 to 90 and pitch angles from -40 to 40 . Both 

yaw and pitch angle were varied at Mach numbers 0. 2 and 0. 5, Yaw angle 
alone was varied at Mach numbers 0. 3, 0. 35, and 0. 4. A description of 
the free jet apparatus is given in Appendix B of Ref. 1. A Prandtl probe 
was used to measure the total pressure at the same radial station as the 
test probe. A conventional strain gauge transducer and scanivalve were 
used to read the pressures from the Prandtl probe, the probe under cali- 
bration and the atmospheric reference. Frequent cross checks to a water 
column manometer were made to insure accuracy. An example of the 
results is shown in Fig. A2, where the pressures measured in yaw surveys 
at five Mach numbers are plotted. 

A. 2 Analysis of Results 

The probe pressure, Pp, at each yaw and pitch angle was reduced 
to a nondime ns ional pressure coefficient Cp defined as 
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— p. - p. 

Cp = 



J/op 



(Al) 



s M 



where Pj = reference impact pressure from the Prandtl probe, Pg = 
static pressure (atmospheric), M = Mach number calculated using refer- 
ence impact pressure and atmospheric static, pres sure, and = 1.4 = 

ratio of specific heats. The values obtained for Cp were plotted against 
the angle of the flow to the probe axis,H^ . The results are shown in 
Fig. A3. For angles that were a combination of pitch and yaw angle, 
the flow angle relative to the probe axis was calculated using the geomet- 
rical relationship. 

= cos ^ (cosO^ cos (3) (A2) 

where 0^ is the yaw angle and 0 is the pitch angle of the flow to the 
cylindrical axis of the probe. Table Al shows the combinations of pitch 
and yaw angles from which the angles in Figures A3g and A3i were 
computed. 

The expression chosen to represent the data at each Mach number 
was = A(sin^ B( - V At least squares algorythm was developed 

for the computer which determined the constant coefficients A, B, 

o 

and N from given data. The algorithm is described in Appendix B. 

Further examination revealed that for this probe the Mach number depend- 

o 

ence could be accounted for by redefining A = A * M* so that from 
Equation Al, the expression 
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Cp = A (sin^ B 



(A3) 



where Cp was defined as 




(A4) 



and A, B, q ^ were constants, described the behavior at all Mach 



numbers tested. This expression was found to hold well as long as the 

yawed probe pressure did not fall below the static pressure, which occurred 

o 

at a flow angle of approximately +60 . Application of this expression (Eq- 
uation A3) in a method for representing the calibration of the Dodge probe 
is described in Appendix D. The expression was used in the present work 
to determine the yaw angle from measurements made with a cylindrical 
impact probe set at different angles to the flow. This application is dis- 
cussed in Appendix B. 

In the pitch angle surveys of the pneumatic probe shown in Fig. A-1 
there was a noticeable assymetry in the distributions thought to be an effect 
caused by the probe stem. The assymetry should be noted in the applica- * 
tion of the third method of calibration for multiple sensor probes described 
in Appendix D. 
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TABLE A-1 ANGLE SETTINGS USED TO COMPUTE FLOW ANGLES IN 

FIGURES A-3g AND A-3i 




FIGURE A-1. PHOTOGRAPH OF 0. 032" DIAMETER IMPACT PROBE 
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FIGURE A-3c. YAW ANGLE SURVEY FOR M = 0. 353 
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FIGURE A-3d. YAW ANGLE SURVEY FOR M = . 399 
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FIGURE A-3e. YAW ANGLE SURVEY FOR M = 0. 500 



79 



FLOW RNKLE 




p^r^N]f^r^r^N!]p^r^r^N3— * 

1 t 1 1 1 1 I i I 1 1 I 1 

PRESSURE CDEFFICIENT 

FIGURE A-3f. PITCH ANGLE SURVEY FOR M = . 204 



80 



FLOW RNGLE 




1 




pressure: coefficient 
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APPENDIX B 



LEAST SQUARES APPROXIMATION OF DATA USING 
Cp = A(sin^B( y - y 

1. Derivation of A, B, N and ^ 

From surveys of a cylindrical impact probe in a known flow, pres- 
sure data is reduced to a set of pressure coefficients C i and flow angles, 

P 

For n sets of points (Cp^, i) the coefficients A,B,N and 



form a least squares fit when the expression 

Q = E (Cpi-A(sin^B( ^ 
i=l 



(Bl) 



is minimized. A minimum value of Q is found when, simultaneously: 

/ \ 



n 



= -2 Z 

i = l 



[Cp.-A(sin^B( V. - y ))^] 

1 o 

[sin^B( y . - yo)]^ 



=0 (B2) 



^ =4AN j; j [Cp.-A(sin^B( y . - . [cosB(y. - y^)]. 



C)B 



i=l 



[(sin^B( yi-y^))^'^]. [sinB(y.-Y^)].[y.-y^] 



n 






1 = 



= 4ANB Z [CppA(sin^B(y.-y^))^]. [cosB(y. -Y^)] 






=0 

(B4) 



[(sin^B(y^-Yo))^"^]. [sinB(y.-y^)] ^ 



<^N = -2A Z 



[Cp.-A(sin^B(y. -Y^))^]. [(sin^B(Y. -y^))H 



i=ll[ln(sin^B(y.-y^))] 



(B5) 



0 (B3) 
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The simultaneous solution of equations B2, B3, B4 md B5 was achieved 
by a numerical procedure which sequentially solved each equation to find a 
value for the respective coefficient while the remaining coefficients were 
held constant* The procedure stepped through the solution of each equation 
in a cyclic manner until the simultaneous solution for A,B,N and 
was found. The solution of each equation for the value of the respective 
coefficient was achieved in an iterative manner using Newton’s method for 
successive approximations. The derivative used to compute the iteration 
interval was defined as the linear slope of the function computed with the 
most recent pair of values for the iterated coefficient. The first two values 
of the function were obtained using a fixed interval to increment the co- 
efficient. The solution was found when the value of the function was less 
than a specified small quantity. 

The simultaneous solution of the four equations was obtained when 
the value of each respective partial derivative computed on the first 
iteration step was less than a specified small quantity. 

No convergence criteria was applied to the sequential solution 
technique. After each individual coefficient was computed its value was 
held constant in the calculations for the other three coefficients. 

The BASIC program ’’ITERATIVE LEAST SQUARES” for this method 
is listed in Table B-1. The computer time required to obtain a simultan- 
eous solution for the coefficients A,B,N and depended on the number 

of data points used and the accuracy of the initial approximation for the 
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coefficients. Experience gained in using this procedure indicated that the 
following guidelines were appropriate: 

a. Critical data points are in the ranges of positive and nega- 
tive angles which define the nearly linear regions of the curve. The three 
nearly linear regions of the curve include, one centered about =0, 

for which C is greater than -0.02, and the two sides with slopes which 
P 

are equal but of opposite sign. The two sides are limited on one extreme 
by the maximum and minimum angles and on the other extreme by the 
angles at which significant curvature begins. Thus six data points, two 
in each of the nearly linear ranges, are required for a reasonable fit. 

b. When selecting additional points consideration should be 
given to the areas where the most accurate characterization of the data 
is required. 

c. The initial estimate of should be the angle, in degrees, 

o 

at the center of the C vsy curve. 

P 

d. The initial approximation of B should be 1.0. 

e. The initial approximation of A should be 0.5 less than the 

minimum C of the data. 

P 

f. The initial value of N should be estimated by dividing the 
angular width, in degrees, of the linear region for which Cp is greater 
than -0.02, by 18 degrees. 
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2 . 



Determination of the Zero Yaw Angle 



Using measurements of Cpj obtained from a cylindrical impact probe 
in a known flow, the coefficients A,B,N and are determined using 

the procedure described above. A,B, and N now characterize the probe 
with respect to Mach number and the actual angle of the flow to the probe 
axis. If is the angle of the probe axis with reference to some mech- 
anical scale, is a measure of the difference between the zero on the 

mechanical scale and the probe angle at which the probe axis is aligned 
with the flow (zero flow angle). This is illustrated in the following 




When A, B, and N are known, measurements taken in an unknown 
flow with the same probe can be used to find which is the angle of 

the flow referenced to the mechanical scale. The pressure insensitivity 
of the cylindrical probes near the zero flow angle and the relationship 
in Equation A2 allows the assumption to be made that for small pitch angles 
the calculated zero flow angle of the probe will be equal to the zero yaw 
angle. 
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In order to determine the zero yaw angl^ measurements of pressure 
are made at a series of probe angles about a reference value. Equal 
magnitudes on either side of the reference are taken (i. e. H^= 0, 

^ ^ i ^4^ max^* procedure to calculate the zero yaw angle 

of the flow is as follows: 

a. Reduce pressure data to values of C . In a known flow 
this is simply to use Equation A4. In an unknown flow the C must be 
approximated since the total and static pressures are not known. A suit- 
able approximation is obtained using the equation 



Cp=A(sin 



N 



P- -P 
1 max 



P +P 

P -f min nmin) 
max ^ ' 



(B6) 



where; A, B and N are previously determined calibration constants, 

P^ = ith impact pressure in the series 

P = maximum P. in the series 
max 1 

Pmin ” minimum P^ in the series 

- next to minimum P in the series 
nmin ^ 

Y = absolute magnitude of the maximum and minimum 

max 

probe angles in the survey. 

This equation scales the pressure data to form a coefficient which has a 
variation similar in shape to the variation of Cp determined in the calibra- 
tion test. 

b. Determine ^ from the data by solving equation B4 for 






while keeping A,B, and N constant. The resulting ^ will be 
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the least squares solution for the zero yaw angle of the flow with respect 
to the mechanical probe angle scale. 

A BASIC program "KAW40" which accepts pressure data and out- 
puts the value of is listed in Table B2. 
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APPENDIX C 



DATA ACQUISITION SOFTWARE 

C-1 BASIC PROGRAM "DATACQ” 

The BASIC program "DATACQ" was written to acquire data using 
the Hewlett-Packard HP 21MX computer and peripherals described in 
Section III of this paper. It is modular in construction and can be used 
by a relatively inexperienced operator to sample real time data in a 
variety of modes. When used interactively as a multipurpose data acquisi- 
tion control program, a question and answer format is used. For special- 
ized data acquisition tasks, it can accept by merging, user subprograms 
which enter the modules (subroutines) in "DATACQ" in a non-interactive 
manner, thus eliminating the need for keyboard entries. 

The key to the flexibility of "DATACQ' is the "Selector Module" 
which provides the control to access all data taking subroutines on command 
A brief "user manual" is provided in this appendix to acquaint the user with 
the operation of "DATACQ". A program listing is provided in Table C-1. 

C-2 USER MANUAL FOR BASIC PROGRAM "DATACQ" 

A limited amount of information is provided in the remarks con- 
tained in the program. This users manual provides instructions for each 
function preprogrammed into "DATACQ". The first step in using "DATACQ' 
is to enter a user program that is written separately then merged into the 
main program. The user program must be written between lines 20 and 
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2000 to avoid overwriting other program lines. The simplest user pro- 
gram that can be used to acquire data is: 

21 GOSUB 5000 

22 STOP 

This program transfers program control to the selector subroutine. 

Subroutine 5000 - Selector Module 

Upon calling subroutine 5000, the query 

ENTER SELECT VARIABLE 
? 

will be displayed. A number corresponding to the desired module in the 
following list should be entered: * 

1 - Initialize calibration coefficient array 

2 - Enter calibration coefficients from console 

3 - On line probe calibration 

4 - Synchronized sampling of data 
6 - Free run sampling of data 

Entry of any number other than the numbers listed here, 0and5, causes the 
query "ENTER SELECT VARIABLE" to be repeated. After a selected 
module has completed its operation "ENTER SELECT VARIABLE" will be 
repeated until a 0 or a 5 is entered, signaling the program execution to 
return to the user program. 

Subroutine 4000 - Initialize Coefficient Array 

This subroutine initializes the calibration coefficients for all data 
channels to have the values Xq = 0 and x^ = 1. As a result, samples of 
data taken on all channels will be displayed directly as voltages. 
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Subroutine 4050 - Enter Calibration Coefficients 



This subroutine allows console entry of the calibration coefficients 
Xq and Xj for all 16 channels. Previous coefficients for the user selected 
channel are lost upon entry of the new coefficients. 

Subroutine 3500 - Time 

If the proper time of day was entered into the computer on startup, 
this subroutine returns with Rl, T2, and T3 set to the time of day in hours, 
minutes past the hour, and seconds past the minute respectively. T is 
set to the time of day in seconds. 

Subroutine 3600 - Time Scheduling 

Caution - This subroutine maybe called only once during "DATACQ” 
execution. A second call will result in an indefinite waiting state in pro- 
gram execution requiring an "abort" command to return computer control 
to the operator. 

The purpose of this subroutine is to allow repetition of a task with 
start time, number of repetitions, and the interval between repetitions 
defined by the user. The task to be repeated must be programmed with 
the first program line number being 100, 300, 500, or 900 and a "Return" 
command as the last program line. Before calling the subroutine, the 
variable V9 must be set to a value equal to the first line number of the 
task to be repeated. 

On completion of the required number of task repetitions, execution 
of a users program will continue with program line 30. 
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Caution - The user must insure that the time interval entered is 



sufficient to complete the task scheduled. Overlapping of scheduled tasks 
is not allowed. 

Subroutine 2000 - RPACE Subroutine 

This subroutine is used for synchronized sampling of a selected data 
channel. Before calling the subroutine the following variables must be 
defined. 

A1 = 32768 + 256 * (Blade Pair) + SAMPLE LOCATION (0-255) 

N1 = number of data samples to be taken at the specified location 

(Maximum value 2000) 

N2 = number of channel to be sampled 

The returned variables and their values are: 

A2 = (15xl0^)/RPM 

C(20, 100) = Data. The first N1 elements of this array (in row 
order) contain the values of the data samples. 

C0 = The numerical average of N1 samples 

Subroutine 2500 - Free run sampling 

This subroutine is used for non- synchronized, "Free Run" sampling 
of a selected channel. Input variables N1 and N2 are defined as in the 
RPACE subroutine. The N1 samples requested are taken at 10 micro- 
second intervals. Output variables C(20, 100) and C0 are as defined in 
the RPACE subroutine. 
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Subroutine 7300 - Subroutine to sample pressure 

This routine will sample channels designated as pressure channels 
and a channel designated as a reference pressure channel in the free run 
mode, and reduce the voltage data on-line to pressures using the expres- 
sions 

Pr ~ ^^o^r (^i)r*(P)r 

and 

Pp = (*o>p ■> (*l)p'(E)p +Pr 

where the subscript, r, refers to the reference channel and subscript p, 
refers to the pressure channel, x^ and x^ are the calibration coefficients 
for the selected channel, E is the ensemble average of the sampled volt- 
ages, is the reference pressure and Pp is the calculated pressure for 
the selected channel. 

Inputs are defined interactively by answering the queries printed 
on the CRT console or, for non-interactive use, values may be pre- 
assigned to the listed variables and the subroutine entered at line 7340. 
The variables which must be defined in the non-interactive mode are: 

E3 -Channel # for the reference pressure or enter 16 for input at 
the console 

E4-Number of free run samples desired on reference pressure 
channel 

E5-Channel number of the pressure channel to be sampled 

E6-Number of free run samples desired of the pressure channel 
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The subroutine prints the pressure, reference pressure, and pres- 
sure channel average voltage on the console CRT. 

C-3 CONTROL PROGRAM FOR COMPRESSOR TEST 

A program was written to be merged into. "DATACQ" which caused 
data to be output on the tape punch and sent to the HP 9830 calculator. The 
program was used for the compressor test described in Section V. A list- 
ing of the control program is given in Table C-2. 
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REM REAL TIME AQUISITION PROGRAM K A WINTERS 1 MAR 78 
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3230 LET N= 
3240 RETURN 
3250 REM 
3260 REM 
3270 REM 







-J 


z 










<a: 


01 










i> 


UJ 










QC 


-J 










UJ 














s 










z 


< 










1 


O) 










Q 


z 










Z 


UJ 










< 


Ul 












3 










UJ 


H» 










s 


UJ 










1 -^ 


X 




o 






H* 






o 








UJ 




o 






H» 


h: 










X 


>— 1 




+ 






C 


M 
















H 






01 


O) 




II 


X 






UJ 




o ^ 






o 


-J 




^ Hi 


H- • 




UJ 






+ 


3 




X 


s: 




^ H 






o 


<c 




H UJ 


01 




01 


to 01 




** -J 


X oi 












oc o 




o 


+ u. 




Hi O 


=3 Z 






o o 




O 


H* O 






o 




H O 


UJ O 




O 


— X 




UJ VO 


cc u 




O 


* UJ 




-J 


01 






CO n 




II 


UJ II 




to 


1 - s 




O A 


Z lO 


^ o 




+ 13 




VO 


I-H H» 


o ^ 


'31 


O z 




O 




* 


ID 


o 




II O 


=3 01 


* 00 


01 


O UJ 




A O 


O UJ 


/-N H» 


o 


o s 




O 


CC H» 


1 


o 


1-4 




O -4 


X :d 


H» O 


1 


* H» 




O * 


:3 z 


1 o 










Ol H-l 


^ o ^ 


CD 






* o 


s 


O O fO 


z 


•^X 




^ o 


II 


o ^ * 


1— 1 


= < 




o o 


Ul CO 


VO fO -• 


-J 


3 Hi 




o o 




fO \ H» 




O 01 JA ^ 






a. ^ 


\H» 1 


s 


Z H ^ 




\ \ 


S 01 


H* ^ 


< 


X H» 




A ^ ^ 


<c cc 


w w 


Ol 


Ul UJ C 7 > • 




H» H Hi 


Ol ID 


M M 


O s 2 




4 > ^ ^ 


O 


z z z 


O O H-i z CX) 


O ^ M Hi 


UJ X 


H - 1 K - 1 


UJ lA H* UJ • VD 


C 30 H Z Z 


S H 


II II II 


s 


: : H* 11 fO 




II H-l 


1 — • 


H» CO »0 Z 


»— 1 


< 7 ^ s-^ 


^o 


-q* 1 1 


H» H» H» H» X 


!-• CQ H» H» H» H-i Z 




H* "T ^ 




UJ =J 


1 ID 


2 Z O O 


O 


Hi Hi 


SSEH'Mt-'J-'SSSECOi-l^-tCLt-'Z 


H 


Luuji-fUiujLtjujujuiujujooscczuja: 


O UJ U. U. 


Qsat-jj-jfxircccciroQ.a.'-i-jH* 


CD -J 1 — I HH 


OOOOOOOOOOOO "0 0000 % 


O lA O O 


o 


cOfO'^iAV^)r*-coc 7 ^oooio^tAr^r*- 


CO CO cn O 


miAtniAintnmiAiAin 




VO VO ^ r-- 


» 0 f 0 i* 0 » 0 i 0 i 0 f 0 i 0 i* 0 i 0 i 0 i 0 i 0 ^ 0 i 0 i 0 i 0 


to to to to 



o 

o 

o 

0 

Os3 

1 

• I 



Ul 



o 

o 



O CQ 

o :3 

O 05 

o o 

<r ^ (S3 

CO + 

o 

II H-l o 

A II 
Oi II 
^ H-I <7^ 

H* r> 



H* 

Cx. LxJ U* 



o m o 

-4 -• OsJ 
r** 

fO lO »o 



O 


o 


o 






o 


o 


o 






^o 


A 


r- 


o 




CQ 


X 


CQ 


to 




ID 


ID 


ID 






Ol 


01 


01 


z 


■q* 


O 


o 


O 


UJ 


Hi 


CD 


CD 


CD 


X 










H 


A 


O 


o 


O 




CO 


O 


o 


O 


o> 


VO 


to 


A 


r^ 


z 


to 


II 


II 


II 


A 




a% 


(7S 


<71 




z X 


3 


Z> 


Z> 


1-4 


o :d 

Z H 


U. 


u. 


u. 


u. 


X UJ 


1-4 


I— I 


1-4 


1-4 


Hi X 




OsJ 


to 




Ol o 


CO 


CO 


CO 


CO 


to <r 




H- 


H» 


r- 


r^ 


to 


to 


to 


to 


rO to 



TABLE C-1 (continued) 



105 



3740 RETURN 



4000 REM CALIBRATION COEF CHANNELS 0-15 ARRAY T(CHN#+l,2) 

4010 REM SUB 4040 INITIALIZES WITH X0=0,X1=1 

4011 REM SUB 4050 INPUTS COEF FROM KEYBOARD 
4040 FOR 1=1 TO 16 
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5190 GOTO 5100 
5200 RETURN 
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7160 GOSUB 2500 
7170 LET X=CO 
7180 LET Y=P6-P7 
7190 PRINT# 6;P6,P7,X 
7200 GOSUB 3050 
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7500 REM SUBROUTINE FOR PACED SAMPLING 

7501 REM AVG VOLTAGE FOR EACH DISCRETE SAMPLE POINT IN ARRAY D(POINT) 

7502 REM INTERACTIVE ENTRY SUB 7500 

7503 REM PREDEFINED ENTRY 7540jJ2=BL PR,J3=SAM SP,J4=1 PER REV,N1,N2 
7510 PRINT "ENTER NUMBER OF SAMPLES, CHAN#" 
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1 




50 GOSUB 5000 
100 REM DATA OUTPUT 

110 PRINT "INPUT ANGLE FOR KULITE PROBE 

111 INPUT A 
IZO PRINT# 4;A 
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TABLE C-2. COMPRESSOR TEST CONTROL PROGRAM 



no 



290 PRINT "DATA TRANSFERED AND PUNCHED 
300 PRINT 
310 GOTO 50 



APPENDIX D 



DERIVATION OF MULTIPLE SENSOR PROBE CALIBRATION 
USING SINGLE IMPACT PROBE CHARACTERISTICS 



Using the characteristics of the cylindrical impact probe discussed 

t 

in Appendix A, which was similar in geometry to the individual sensors of 
the Dodge probe, a calibration was derived for the Dodge probe. 

The following definitions were used: 



1 



23 



4 

H^i 

a 

0 

M 



- Pressure at the Dodge Pj^ sensor 

- Average of the Dodge P 2 s^nd P 2 sensor pressures 

- Pressure at the Dodge P^ sensor 

- The angle of the flow with respect to the sensor axis 

- Yaw angle with respect to axial direction 

- Pitch angle with respect to the plane normal to the probe 
shaft 

- Mach number 



Pg - Static pressure 

P^ - Total pressure 

^ - Ratio of specific heats (^=1.4 was used) 

In Appendix A it was shown that, for a single impact tube, the 
pressure coefficient defined as 



C = 
P 



Pj-Pt 
T PgM' 



( 1 ) 



was closely approximated by the expression 
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! 



( 2 ) 



Cp = A [sin^B{ y i- y^)]^ 
where the flow angle was determined from pitch and yaw angle by 
y .=cos"^{cos Ct JCOS 0j) (3) 

For the Dodge probe, when rotated to balance the pressures at sensors 
and P^, the following table lists the values of the angles for each of 
the sensors 



a. 

1 


Dodge 

1 


sensor 

2-3 
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0 


a 


0 


0i 


0 


0 


0 -^ 




0 cos 


~^(COs(y COS0) 


0-0 



where 

0^ = angle between axis of sensor 2 or 3 and sensor 1 

0 = angle between axes of sensor 4 and sensor 1 



The quantities which have been used previously to represent the probe 

1 

characteristics are: 





-P 



-P 



4 



23 



and 




P 




Using equation (2) in equation (4) 



(4) 



( 5 ) 



J , °pl ~ °p4 
^pl“^P23 



[sin^B yi]^ - [sin^B y 4 ] ^ - 
[sin^B y - [sin^B H^23^^ 
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Using equation (3), 



[sin^B0]^ - [sin^B(0-0)]^ 

Q = 



[sin^B0]^”[sin^B(cos " ^^cos CX cos 0^ )] 



N 



( 6 ) 



Note that if the pressure coefficient, Cp is independent of Mach number, 
then equation (6) gives pitch angle (0) in terms of sensor measurements ( ^ ). 
There must be no interference between sensors, however for equation 6 
to describe the probe characteristics correctly. Using Equation (5) 

^ ^ ^P . 23 _ A[sin^B0]^ - A [sin^B(cos’^ { cos cos 0) )]^ 

P _ -D.. .r . 7 „wiN 5 I 'V 1 (7) 









A[sin2Ba]'^ - -J— , [1 



7 



M' 



Defining 



X = 



_ 2 [ 1 + j.r } . ^ 

2 2 



7 M 

equation (7) may be written 



( 8 ) 



X = A[(-^ -1) (sin^ (b 0) )^- A (sin^B(cos*^ {cosC^! cos 0> ))^] 

P J2 

If the pitch angle is known, equation (9) gives the Mach number, M, for a 
measured value of , 

Alternatively if we define 



(9) 



6-0- 



^pl ~ ^p4 

Cpl ^ ^ 



( 10 ) 



then equation (9) becomes; 



- A [(-^ -1) (sin (b 0) ) _ ^ )^j 



113 



Equations (6), (8), (9) and (11) were used to derive a calibration 
for the Dodge probe. The coefficients A, B, and N which represent the 
characteristics of the single impact tube similar to the sensors of the 
Dodge probe were determined as in Appendix B. 

From measurements taken in a known flow at zero pitch angle 
(0 = 0) with the Dodge probe the quantities J3 , 'Y and (5 were calculated 
using equations 4, 5 and 10. 

The angle C( was computed from equation 9. 

First, X was calculated from equation (8) using the measured 
Mach numbers. Then from equation (9), 

^ - ^(sin^ BCX ) 



or 



sin 

a = 



0 



2N ] 



B 



The angle ^ was computed similarly from equation 11: 

X 6 , — 1 — T 

"lU — Z-K ) 2N ] 



sin 



B 



Having established 0^ and ^ from measurements the procedure to com- 
pute velocity using values of ^ , 'Y ^.nd (5 measured in an un- 

known flow was: 



1 ) 

2 ) 

3) 



Assume 0=0 



- / 



Calculate the right hand side of equation (> - Y 



If 



7 - 7 ' 



~ 0 then a solution for 0 was found. Proceed 



to calculate M. 
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I 




4) Calculate 



d7ld0 

7 - 7 ' 

5) 

6) Return to step 2 with the new value of 0 

Repeat the iteration .until a solution for 0 was found 
To calculate M: 

7) Calculate the right hand sides of the equations (9) and (11). 

Compute the average = x . 

8) Assume a starting value for M 

9) Calculate x using equation (8) 

10) If |k “ 3<j ~ 0 then a solution was found and the 

unknown flow was determined 

11) Calculate cf x/ dM 

X — )T 

12) Set M = M + ^ x/d M 

13) Return to step 9 with the new value of M 

Repeat the iteration until a solution for M was found. 

This method (method III) for calibration of the Dodge probe was tested using 
data previously obtained with the Dodge probe in a known flow. Calibration 
methods I and II were used to reduce the same data to pitch angle and 
velocity magnitude, (Each data point corresponded to values of pitch 
angle and velocity magnitude which were controlled and known during the 
calibration test.) Table D-1 is a comparison of the accuracies achieved 
using the three methods. 
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It can be seen that method III was the least accurate method of 



representing the calibration of the Dodge probe. However, the mutual 
interference of the probe tips is probably the cause of the large error. In 
the 2 probe system of velocity measurement there can be no interference 
between probes and hence this method should be further considered for 
that application. 
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MACH NO. PITCH ANGLE 



CALIBRATION 


AVG. 


MAX 


AVG. 


MAX 


MACH # 


METHOD 


ERR % 


ERR % 


ERR DEG 


ERR DEG 


. 35 


I 


2. 6 


3. 2 


0. 5 


-P. 9 


.35 


II 


0. 5 


-1. 2 


0. 6 


-1. 0 


.35 


ni 


4. 5 


7.3 


0. 8 


-1. 1 


.44 


I 


1. 5 


-2. 2 


0. 4 


-1. 0 


.44 


II 


0. 3 


-0. 6 


0.4 


-0. 6 


.44 


m 


4. 7 


8. 3 


1.0 


-1.4 


.49 


I 


2. 0 


-4. 2 


0. 6 


-1. 8 


.49 


II 


0. 5 


0. 9 


0.4 


0. 7 


.49 


III 


4. 9 


8.4 


1.4 


-1. 6 



Average and maximum values of the errors are shown for the five pitch 
angle settings -5°, 0°, 5°, 10° cS< 15 . 



TABLE D-1 SUMMARY OF ERRORS IN THE USE OF 
THREE CALIBRATION METHODS FOR THE DODGE PROBE 
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I 



I 



I 



1 

I 

i 
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